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Abstract 



The nuclear option at THERA provides an ideal and unique opportunity 



^ ■ to investigate the black body limit (BBL) of high energy Deep Inelastic Scat- 

S-T 



tering of highly virtual photons off heavy nuclear targets and thereby probe 
QCD in a new regime. At high enough energies, whichever hadronic configu- 
ration the photon fluctuates into, the interaction at a given impact parameter 
with the heavy nuclear target will eventually reach its geometrical limit cor- 
responding to the scattering from a black disk. An attractive feature of the 
BBL regime for a large nucleus is that the interaction is strong although the 
coupling constant is small. Predictions for longitudinal and transverse dis- 
tributions of the leading hadrons in inclusive and diffractive channels and 
exclusive vector meson production are found to be strikingly different in BBL 
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and the leading twist approximation. In particular, the multiplicity of leading 
hadrons in the current fragmentation region is strongly suppressed, while the 
cross section of diffractive vector meson production by longitudinally polar- 
ized photons is oc 1/Q 2 - The x and Q 2 ranges where BBL maybe reached are 
calculated for the interaction of color triplet and color octet dipoles at central 
impact parameters. We conclude that for heavy nuclei THERA would reach 
deep into BBL for a wide range of Q 2 . The connection between hard diffrac- 
tion and leading twist nuclear shadowing is analysed. It is demonstrated that 
the current leading twist analyses of the HERA diffractive data lead to similar 
and very large leading twist shadowing for gluons. Several processes sensitive 
to the amount of nuclear shadowing for gluons are discussed. 



I. EXECUTIVE SUMMARY 

The nuclear option at THERA provides an ideal and unique opportunity to investigate 
the black body limit (BBL) of high energy Deep Inelastic Scattering (DIS) of highly virtual 
photons off heavy nuclear targets and thereby probe Quantum Chromodynamics (QCD) 
in a new regime. At high enough energies, whichever hadronic configuration the photon 
fluctuates into, the interaction at given impact parameter with the heavy nuclear target 
will eventually reach its geometrical limit corresponding to the scattering from a black disk. 
An attractive likely feature of the approach to the BBL regime for a large nucleus is that 
the interaction is strong although the coupling constant is small. Predictions for many 
observables are qualitatively different within the BBL from those of the standard leading 
twist DGLAP |l[ regime. The most striking predictions for inclusive observables include: 
(i) the nuclear structure functions behave as F^ oc ttR\Q 2 \yi{1 / Am nRa%), with a 1/x 
dependence slower than that predicted by DGLAP (this is in contrast with the nucleon case 
in which oc Q 2 kR 2 n (1 + c^ln 2 1/x) \nl/x whose x dependence is comparable with the 
DGLAP prediction, here Ra,Rn are radii of nucleus and nucleon), (ii) the ratio of nuclear 
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and nucleoli structure functions decreases with 1/x, F^-jF^ oc 1/(1 + In 2 1/x). Major 
new features for the final states in DIS are: (i) a much softer leading hadron spectra with 
enhanced production of high p t jets in inclusive and diffractive processes, (ii) the inclusive 
diffractive cross section reaches nearly half of the total cross section, (iii) a weakly Q 2 - 
dependent parameter-free cross section for exclusive vector meson electroproduction. 

It is possible that the black limit will be reached for certain gluon-induced DIS processes 
off nucleons within THERA kinematics. However, nuclei have two very clear and useful 
advantages over nucleon targets. Firstly, because the target is extremely Lorentz contracted, 
the partons of individual nucleons are piled up on one another. This means that for a 
given photon-target centre of mass energy the black limit will be reached much earlier in 
nuclei. Secondly, scattering from the edge of nucleons, which remains grey even at very high 
energies, (and competes with the BBL contribution) is heavily suppressed for large nuclei. 
The purpose of the rest of this report is to expand on the details for the interested reader. 

II. INTRODUCTION 

The standard theoretical approach to inclusive Deep Inelastic Scattering (DIS) of a 
virtual-photon off a target, j*(q) + T(p) —>■ X, uses a factorization theorem |2j to prove 
that in the limit of large photon virtuality (q 2 = —Q 2 <C 0), leading twist (LT) terms domi- 
nate in the Operator Product Expansion || of the relevant hadronic matrix elements. This 
corresponds to a factorization of short and long distance contributions and leads to an at- 
tractive physical picture of the virtual photon scattering off electrically charged "point-like" 
partonic constituents within the target, with its associated approximate Bjorken [|J scaling 
(which states that structure functions, for example, F^ and F 2 , only depend on the ratio 
x = Q 2 /2p ■ q). The application of the renormalization group to perturbative Quantum 
Chromodynamics (pQCD), embodied in the DGLAP |l| evolution equations, corrects this 
naive picture and leads to logarithmic scaling violations. These equations describe how the 
long-distance boundary conditions, the parton distributions of the quarks and gluons, (and 



3 



ultimately Fl and F2) change logarithmically with Q 2 . 

It is now understood that this decomposition in terms of twists should become inappli- 
cable in pQCD at sufficiently small x (i.e. the higher twist (HT) terms, which are formally 
suppressed by positive integer powers of 1/Q 2 , become numerically important). However, 
in practice, it is very difficult to distinguish HT effects in the inclusive structure functions 
from uncertainties in the input parton distributions for a nucleon target, because DIS off the 
nucleon edge masks the importance of the HT effects. We will explain that the theoretical 
description of DIS off heavy nuclei is significantly simpler than that for a nucleon target. 
The energy and Q 2 dependence of structure functions of heavy nuclei, and specific properties 
of final states, can serve as clear signatures of a new regime, the black body limit (BBL) 
regime, where significant contributions to the dynamics of the strong interactions from all 
twists are expected. 

For the past decade, pQCD has been used extremely successfully to describe inclusive 
DIS processes (see e.g. |J). Many new phenomena specific to small-x physics have been 
discovered experimentally and explained by pQCD. New QCD factorization theorems for 
certain hard processes (such as diffractive scattering ||, exclusive vector meson production 
|7].|S|] and Deeply Virtual Compton Scattering 0) have been proven. These factorization 
theorems were used to predict and describe qualitatively new small-x hard scattering phe- 
nomena. In particular, hard exclusive processes for both nucleon and nuclear targets, with 
properties strikingly different from those familiar from soft hadronic processes, were pre- 
dicted and observed (several experimental |10[ and theoretical [|lT|-|i"3|| reviews are available 



in the literature). These include: (i) the complete transparency of nuclear matter in coher- 
ent pion diffraction into two jets off a nuclear target (predicted in |14| and later observed 
in ||15|| ) and in exclusive photoproduction of J/ip mesons off nuclei |16[], (ii) a fast increase 



with increasing energy of cross sections of hard exclusive processes ||17|| : high Q 2 exclusive 
p, (f), J I ip meson production |7|JT8|-pT| and exclusive photoproduction of J / ip and T mesons 
P5| and low-x exclusive photon production in DIS (Deeply Virtual Compton Scattering 
), (iii) an almost universal and small slope of the t-dependence of hard exclusive processes 



(predicted in [p5|D, (iv) a small and process/scale dependent rate of change of the slope of 
t-dependence of exclusive processes with energy, a', (v) confirmation that QCD is flavour 
blind in hard (high Q 2 ) exclusive vector meson (p and mesons) electroproduction. These 
theoretical and experimental discoveries have moved the focus of investigations in QCD to 
new frontiers. 

A seminal and extremely important experimental observation is the fast increase with 
1/x of the proton structure function, F 2 , at small x discovered by the HI p6| and ZEUS 



27j collaborations at HERA. The implied fast increase of the proton parton densities raises 



several challenging questions for DIS at even higher energies: 

• Will it remain sufficient to use the leading twist DGLAP evolution equation (to a given 
accuracy in logarithms in Q 2 ) or will it become necessary to resum p8|-p0| the large 
logs in energy that appear at each order in the perturbative expansion (i.e. for the 
regime in which a s (Q 2 ) ln(a;o/a;) ~ 0(1)) ? 

• Does diffusion in parton transverse momenta lead to a breakdown of (collinear) fac- 
torization for leading twist pQCD, and a non-trivial mixing of perturbative and non- 
perturbative effects, in a kinematic regime in Q 2 which is assumed to be safe at higher 
x ? 

• Will the increase of the parton distributions lead to an experimentally-accessible new 
pQCD regime, where the coupling constant is small but the interaction is strong ? If 
so, how can one distinguish between this new regime and the standard one in which 
the LT DGLAP equations are applicable ? 

• Will the structure functions and cross sections of hard exclusive processes continue 
to increase with energy or will their growth be tamed to avoid violation of unitarity 
of the S-matrix (applied to the hard interactions of the hadronic fluctuations of the 
photon) ? 
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• Since the same QCD factorization theorems which lead to successful description of hard 
processes predict a rapid increase of HT effects with increasing energy, how important 
will the latter effects be for the interpretation of the small- a; data ? 

In the following we explain that one has a good chance to answer these questions by study- 
ing electron-nucleus collisions at THERA. We will formulate several theoretical predictions, 
based mostly on general properties of space-time evolution of high-energy scattering, which 
can serve as signals of new hard QCD phenomena. 

DIS at high energies is most easily formulated in the target rest frame (this approach 
has become known as the dipole picture for reasons that will become apparent below). In 
this frame, the photon fluctuates into a hadronic system (a quark anti-quark colour triplet 
dipole, to lowest order in a s ) a long distance, l/2m T x 3> Rt, upstream of the target. This 
system interacts briefly with the target and eventually the hadronic final state is formed. 
The formation time of the initial state hadronic fluctuations, and of the final state, is typ- 
ically much longer than the interaction time with the target. This allows the process to 
be factorized into a wavefunction for the formation of the initial state, an interaction cross 
section with the target, and a wavefunction describing the formation of the final state. In 
this frame, the standard LT result of pQCD, at leading log accuracy, may be re-expressed 



using the following well-known interaction cross-section |p^,|18|j!9|,pl| for the scattering of 
the qq dipole with the target (for an explicit derivation see 



7T 2 

OpQCD^La;) = — d 2 ± a s (Q 2 ) xG T (x,Q 2 ) . (1) 

which is a generalization of the two-gluon exchange model of |33| (see also references in |fL2|j ). 
Here d± is the transverse diameter of the dipole, xG T is the leading-log gluon distribution 
of the target sampled at a scale Q 2 = X/d 2 ± , with A a logarithmic function of Q 2 . There are 
of course corrections, involving more complicated partonic configurations such as qqg, which 
in general can not be described (except in special circumstances) using the interactions of 
dipoles of a given colour configuration. However, the description of small-x processes in 
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terms of an interaction cross section for the hadronic fluctuations of the photon remains a 
useful one. 

For heavy nuclear targets, it is legitimate to neglect a logarithmic increase of the effective 
radius of the target with increasing energy. In this case there is, of course, a natural geomet- 
rical upper limit to the size of this inelastic interaction cross section given approximately by 
the transverse cross-sectional area of the target: 

OpQCD < Oblack= ttRt . (2) 

It seems plausible that small-x physics of the next generation of accelerators will be a 
combination of two generic phenomena, generalised colour transparency and complete opacity 
(absorption), which correspond to two different kinematical limits. In one limit, when x is 
fixed and Q 2 — ■> oo, hard processes exhibit numerous colour transparency phenomena, which 
can be described in terms of various generalisations of the QCD factorization theorem. In 
another limit, when Q 2 is fixed and x — > 0, the interaction cross section increases with 
decreasing x (increasing energy), which leads to increased absorption by the target. Hence, 
in this kinematical limit, the approximation of complete absorption, which could naturally 
be termed "the black body limit (BBL)", seems to be a promising guide to predict new 
distinctive hard QCD phenomena. Colour Transparency is a generic name used to describe 
the fact that small colourless configurations do not interact strongly with hadrons, which is 
reflected by the inclusion of the transverse diameter squared, d\, in (JJ). However, for a fixed 
dipole size (and neglecting for a moment an increase with increasing energy of the radius 
of a target which can be justified for heavy nuclear targets only), if the gluon density of 
the target continues to rise with increasing 1/x, the interaction cross section will eventually 
reach its geometrical limit (cf. (|2])). This is basically what is meant by complete opacity. 

The validity of the BBL can be justified for DIS processes off heavy nuclear targets since 
contributions of colour transparency, and various peripheral phenomena (which complicate 
studies of BBL for nucleon targets) are suppressed. In the BBL, the nuclear structure 
function F^(x, Q 2 ) is predicted to behave as R\Q 2 In 1/x with decreasing x, which is slower 
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than the x-behaviour predicted by the DGLAP evolution equations. Since the nucleon 
is "gray" rather than "black" as seen by the incoming high-energy fluctuation, peripheral 
effects are important and contribute to the total electron-nucleon DIS cross section (these 
effects lead to an increase of the essential impact factors with increasing energy). The nucleon 
structure function (gluon distribution) behave as (xG N ) oc Q 2 R 2 N In 3 1/x in the black 
body limit. In both the nuclear and nucleon cases, the expression for the total cross sections 
is a direct generalization of the Froissart unitarity limit to DIS (see subsection |ITl|) . One 
has to distinguish two features of these results. Firstly, the predicted Q 2 -independence of 
total cross sections grossly violates Bjorken scaling and serves as a very clean signal of the 
onset of the BBL regime. Secondly, the predicted energy (or x) dependence of the total 
electron-nucleon (but not electron-nucleus) cross section (nucleon structure functions ) can 
probably be accommodated within the DGLAP LT pQCD framework by a suitable subtle 
re-tuning of the initial conditions (this has been the pattern in the era of HERA). Hence, 
it would be difficult to distinguish between the DGLAP and black body approximations 
and predictions thereof for the nucleon structure functions. A much more promising way to 
identify the BBL regime is to investigate the difference between the x dependence of nucleon 
structure function and the structure functions of heavy nuclei, as well as final states in DIS 
on nucleons and nuclei for which the BBL predicts several rather striking new phenomena. 

By fixing Q 2 and decreasing x, one proceeds from the DGLAP regime to the BBL regime. 
Theoretical models attempting to quantify deviations from the DGLAP approximation and 
define the kinematical boundaries of the applicability of the DGLAP equations agree that 
the effects, which are responsible for the inapplicability of the DGLAP dynamics, are en- 
hanced for nuclear targets by approximately a factor of A 1 / 3 if the approximation is made 



that the HT effects and nuclear shadowing of parton distributions are a correction p5 
Neglecting non-perturbative nuclear shadowing of the gluon distribution, one can estimate 
the ratio of the critical x at which the BBL regime sets in for DIS on nuclei and nucleon, 
as xa(bbl)/xat(bbl) « (A/2) 1//( ^ 3ri ^ (the factor n comes from the energy dependence of the 
nucleon structure function F% in the form F<^ oc x~ n ). The factor 1/2 in this relation ac- 
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counts for the difference between the nucleon radius, = 0.85 fm, and a mean inter-nucleon 
distance in nuclei, r NN = 1.7 fm. 

Note that the complication of taking non-perturbative nuclear shadowing in nuclear 
parton distributions into account slightly reduces the advantage of using nuclear targets 
over nucleon ones for studying higher twist effects. 

Since the DGLAP QCD evolution equations are leading twist equations, any violations 
of the DGLAP equation would signal a non-negligible role of higher twist effects. Therefore, 
a unique feature of electron-nucleus collisions is that by studying DIS on nuclear targets one 
can not only amplify the higher twist effects but also study them as a function of the target 
thickness (by varying nucleon number A). 

In order to summarize the above discussion, we would like to re-emphasize that using 
nuclear beams in THERA has significant advantages over using a proton beam in studying 
new QCD phenomena such as the perturbative regime of very high parton densities and 
the transition from the DGLAP regime to black body scattering^ Since these phenomena 
depend on the nucleon number, A, measurements of various observables as a function of A 
will be necessary. Fortunately this will be practical at THERA. The BBL interpretation of 
many phenomena, including the difference between the x-dependence of heavy nucleus and 
nucleon structure functions and some properties of final states, does not require knowledge of 
nuclear effects. It is clear that a consistent and complete interpretation of some experimental 
data would require reliable information about traditional nuclear effects at small x, primarily 
nuclear shadowing. However, as we shall explain in detail in Sect. [V], using the profound 
connection between nuclear shadowing in inclusive electron-nucleus DIS and hard electron- 
proton diffraction, the QCD factorization theorem and modern fits to the diffractive data, 

1 Formally, the notion of parton densities is defined in the leading twist approximation only. Thus, 
when discussing physics beyond the DGLAP approximation, one should be cautious and explicit 
in using the term "parton density" . 
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one can significantly reduce uncertainties in the predictions of the leading twist part of 
nuclear parton densities. 

The rest of this section is organized as follows. The black body limit is discussed in 
Sect. [TI1|. We demonstrate that predictions made in the BBL are strikingly different from 
those obtained within the DGLAP approximation, especially for final states in DIS diffrac- 
tion on nuclei. Precursors of the BBL can be studied by analysing the departure from the 
leading twist DGLAP evolution equation using the unitarity of the S'-matrix for hard inter- 
actions as a guide. In Sect. |TV| , we discuss the relevant effects for a range of nuclei at central 
impact parameters. It is found that our predictions are sensitive to the amount of leading 
twist nuclear shadowing, which is considered in Sect. |V|. It is also demonstrated that in a 
wide range of x and Q 2 , leading twist shadowing dominates over higher twist shadowing, 
used frequently in eikonal-type models. After briefly considering experimental requirements 
in Sect. [VI|, we conclude in Sect. [VII . 



III. THE BLACK BODY LIMIT AND ITS SIGNALS IN DIS FINAL STATES 

The leading twist (LT) approximation of perturbative QCD successfully describes DIS 
of photons on hadronic targets || and predicts a rapid increase Q of structure functions with 
increasing energy (decreasing x). Within the LT approximation, there is no mechanism 
which would slow down or tame the rapid growth of the structure functions 0. 

However, it is clear that a rapid power-like growth of the structure functions (at a given 
impact parameter) cannot continue forever [] and, hence, the LT approximation must be 



2 This can be either introduced by hand at the initial evolution scale or generated by QCD evolution 



3 One can show by direct pQCD calculations that the taming (shadowing) of parton densities at 
small x considered in eikonal and parton recombination models is a higher twist effect. 
4 At the same time, conventional nucleon structure functions integrated over all impact parameters 
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violated at very small x. Therefore, some higher-twist mechanism is required to explain the 
taming of the structure functions. A practical and almost model-independent approach to 
the taming of the structure function is based on the analysis of the unitarity of the S'-matrix 
for the interaction of spatially small quark-gluon wave packets. 

A. Unitarity and the black body limit 

The most direct way to understand the constraints which are imposed by unitarity of 
the scattering matrix (for hadronic fluctuations of the virtual photon) on the nucleon and 
nucleus structure functions is to use the impact-parameter representation for the scattering 
amplitude. As explained in the introduction, in the target rest frame the incoming photon 
interacts with the target via its partonic fluctuations. Since the interaction time is much 
shorter than the life-time of the fluctuations at small x, the DIS amplitude can be factorised 
in three factors: one describing the formation of the fluctuation, another - the hard interac- 
tion with the target, and a final factor describing the formation of the hadronic final state. 
The qq dipole is the dominant fluctuation at short transverse distances, and we consider its 
interaction primarily in what follows. 

Within this high-energy factorization framework the structure functions, F L and F T , may 
be written in the following simple way 

F L , T (x,Q 2 ) = [ dzd 2 d ± \r L , T ^dl,Q 2 )\ 2 lmA{S ' t = ^ dl) , (3) 
J s 

where iP1t( z , d\, Q 2 ) are the light-cone wavefunctions of the longitudinally and transversely 

polarised virtual photon, respectively, z is the photon momentum fraction carried by one of 

the dipole constituents, s is the invariant energy of the dipole-target system (s — (P + q) 2 

for the qq dipole), t ~ — \q'±\ 2 is the momentum transfer and d± is the dipole's transverse 

diameter. 



may increase with increasing energy even faster (according to the oc In' 1/x law) than predicted 
by DGLAP. 
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The dipole-target scattering amplitude A(s, t, djj can be expressed via the corresponding 
amplitude, f(s,b,d 2 ± ), in impact-parameter, b, space {b = b± is Fourier conjugate to q' ± ): 

A{s, t, r 2 ) =2sJ d 2 b e iq '*f{s, b, r 2 ) . (4) 

The scattering amplitude f(s, b, d 2 ± ) is related to the total cross section for the scattering of 
dipoles of fixed transverse diameter, d±, by the optical theorem: 

crtotMi) = 2 J d 2 blmf(s,b,d 2 ± ) . (5) 

The unitarity of the scattering S-matrix |37]] (S a b = 8 a b + ifab, Sl c S c b = 5 a b) imposes the 
following conditions on the scattering amplitudes (for each b, a continuum analogy of 
angular momentum in a partial wave analysis in non-relativistic quantum mechanics): 

Imf aa ( S , b, d\) = ^(fljaa + K,c^afLfca) , (6) 

where "a,c" are labels for definite states (on the right hand side we suppress the variables 
s, b, d±_ for clarity). The first and second terms on the right hand side of @ involve elastic 
and inelastic final states, respectively. 

The black body limit (BBL) assumes that: (i) configurations with the impact parameters 
satisfying b 2 < ^ ax are completely absorbed by target, i.e. the elastic matrix elements 
of the S'-matrix are zero for those impact parameters: S aa (b < b m(iX ) = 5 aa + if aa (b < 
bm-Ax) = 0. This implies that Imf aa (s,b < max , d\) = 1, and leads to \f aa (s, b, d 2 L )\ 2 = 
S Cj c^a|/ac('S, b, <ij_)| 2 that is to the equality of the elastic and inelastic contributions to the 
total cross section, (ii) the region b 2 < &^ ax gives the dominant contribution to the scattering 

5 A rather straightforward way to deduce unitarity of the S'-matrix is to consider amplitudes for 
the scattering of bound states of fictitious heavy quarks Q and the prove that, for sufficiently large 
mass of the quarks Mq and small x, dipole scattering would dominate. In order to complete the 
derivation, one also needs to use the fact that QCD is flavour blind if the resolution scales are 
chosen appropriately. 
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amplitude (|j). The great advantage of the BBL is that calculations of the amplitudes of some 
small-a; processes do not require specific model assumptions. Moreover, as will be discussed 
later, the BBL approximation seems to be realistic for DIS on heavy nuclear targets. 

The dependence of f(s, b, on the impact parameter, b, at large b follows from analytic 
properties of the scattering amplitude A(s,t,d 2 ± ) in t-plane. A simple analysis leads to 
f(s, b, d 2 ± ) = c exp(—fj,b) at large b, where c oc xG T (x, d 2 ^) oc l/x n . Hence, following Froissart 
|34] , |3~7H we may evaluate the maximal impact parameter characterising the black body limit. 
One obtains fe max oc l//i 2 ln 2 l/x for a nucleon target, and 6 max = R\ for a heavy nuclear 
target (Ra being the radius of the nucleus). Note that the difference between 6 max for the 
nucleon and nuclear case reflects the fact that the nucleon target is not a homogeneous 
sphere but rather an object with an extended diffuse edge. 

Using these relationships, for the total dipole-nucleon scattering cross section (see also 
(H)) in the BBL approximation, one obtains 

*tot(s, d\) = 2tt(R 2 n + Ac N In 2 l/x) . (7) 

In addition to the total cross section, one can examine the t-dependence of the cross section, 
corresponding to the amplitude in (f|), defined through its slope, B: 

dlna(s,t) o,«/ lTl1/ fd 2 bb 2 f(s,b,dj) 

— dT~ ^ = B = B « + 1/X = 2Jd%f(s,b,dl) ' (8) 

where a' ^ = d (-B/2) /d In l/x. We use above that the t-dependence of the real and imaginary 

parts of the amplitude is the same. In the BBL, one obtains 

5=^ = ^ + c r ln 2 l/:r, 
4 4 

a cS = ct hi l/x . (9) 

Here Rt is the radius of the target, ct is a factor which is similar for nucleon and nucleus 
targets. Hence for practical purposes, one can neglect ct for heavy nuclear targets. 
Finally, the proton structure function (at fixed Q 2 ) in the BBL reads 

o E>2 a 2 

F%(x, Q 2 ) oc £ e 2 Q 2 ^f (1 + || In 2 x /x) In l/x . (10) 
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Here the sum is taken over electric charges of active quark flavours i. Note that the 
additional factor of ml/a; in fllOD as compared to (|7|) is due to the contribution of large 
masses resulting from the singular nature of the photon wavefunction. This reflects the 
logarithmic divergence of renormalization coupling constant for the electric charge (cf. ( [3~3|) 
below). As one sees from (0), general principles of QCD and, in particular, the BBL 
approximation, do not exclude a fast (In 3 1 jx) increase of the structure functions of a nucleon 
at x — > 0. In addition, the contribution from dipoles with the impact parameters larger 
than femax or sufficiently small dipoles (which are assumed to give a small contribution in 
the BBL approximation) should continue to increase with increasing energy as dictated by 
the DGLAP approximation of QCD. Thus, in practice, it would be very difficult, or even 
impossible, to distinguish the BBL prediction (|I]) from a similarly rapid growth predicted 
by the DGLAP equation and to search for saturation effects in inclusive nucleon structure 
functions. Hence, one should turn to DIS on nuclear targets in order to search for distinct 
signals of BBL dynamics. 

The use of nuclei has two clear advantages. Firstly, scattering at large impact parameters, 
where the interaction is far from the BBL over a wide range of energies (the edge effects), is 
suppressed by the factor R n /Ra- Secondly, in a broad range of impact parameters, b < Ra, 
the nuclear thickness is practically 6-independent and much larger than in a nucleon. Hence, 
DIS on sufficiently heavy nuclei can serve as a good testing ground for the application of the 
BBL and will allow us to pinpoint some distinctive features of it. For example, as follows 
from the above discussion, the unitarity of S"-matrix significantly tames the rapid grows of 
the nuclear structure functions and predicts the unitarity limit 

F 2 A (x, Q 2 ) = £ e 2 Q 2 ^ In l/4m N R A x . (11) 

Over the last few years a number of models, using the infinite momentum frame, were 
suggested in order to explain the dynamics of DIS at small x by building the nuclear wave 
function from large gluon fields and assuming a certain saturation of the parton densities 



(for the review and references see e.g., [pq| ). In many respects, these models and the BBL 
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approximation are similar. 



The BBL in DIS from a heavy nucleus at small x was first considered by Gribov [|39 
before the discovery of QCD. Gribov assumed that each hadronic fluctuation of the vir- 
tual photon interacts with the target nucleus with the same maximal strength allowed by 
unitarity. Such an assumption, supported by the observed cross sections of hadron-nucleus 
interactions, was natural for understanding the dynamics of soft strong interactions in mod- 
els predating QCD. 

Thus, in the black body limit, DIS on nuclei is dominated by the dissociation of the 
incoming virtual photon into its hadronic fluctuations, which subsequently interact with the 
target with the same scattering cross sections 2ttR 2 a , and then hadronize into final states with 
mass M. The transverse and longitudinal nuclear structure functions may be conveniently 
formulated as an integral over produced masses 

min 



Im* . 12vr 3 (M 2 + Q 2 ) 2 ' 

mm v ' 

F A (x O 2 ) - C / M ™ X dM 2 Q 4 P( M2 ) 

b L [x, Q ) - L dM ^ , (13) 

mm v ' 

where p(M 2 ) = a e+e ~^ hadrons (M 2 )/a e+e '^ + ^ (M 2 ). In the BBL, the coefficient C = 1. 
The upper cutoff, M max <C W 2 ~ 2q m N , comes from the nuclear form factor: 

^min-^4 _ + Q 2 ) 2 D 2 /o ^ 2 



R 2 A /3 w m 2 N x 2 ^ < 1 . (14) 



3 Aq 2 3 

The key element of the derivation of ( |T3| ) is the observation that in the BBL, as a result 
of orthogonality of the wave functions of the eigenstates of QCD Hamiltonian with differ- 
ent energies, the non-diagonal transitions between states with different M 2 are zero. This 
enables one to write the structure functions as a single dispersive integral as is done in 
(|T3|). Since flT3| ) leads to a cross section for 7*A scattering grossly violating Bjorken scaling 
(er^ot (x, Q 2 ) oc ln(l/a;) instead of oc 1/Q 2 ), the BBL has been considered for some time 
to be an artifact of the pre-QCD physics. This is especially so since, within the parton 
model, the aligned jet model removed this gross scaling violation |40[| . However, for con- 



sistency of the target rest frame and infinite momentum frame descriptions, an exponential 
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suppression with decreasing d\ of the cross section of interaction of small-size configura- 
tions with hadrons was assumed to be required. This fact was subsequently explained and 
understood in terms of the QCD factorization theorem for the scattering of small dipoles, 
colour neutrality of the dipole and asymptotic freedom for hard processes in QCD (and the 
"suppression" was realised only to be only a single power in d 2 ± , cf. ([[])). 

In perturbative QCD, the dipole-target cross section ([!]), rapidly increases with increasing 
energy since the gluon density rapidly increases with decreasing x. Hence, if the increase of 
the interaction cross section is not tamed by some mechanism, it will reach values expected 
for the black body limit (i.e. tens of mb, cf. (§)). 

Properties of the BBL in QCD are somewhat different from those within the Gribov 
picture due to a significant probability of smaller than average size configurations in the 
photon wave function, for which the conditions of the black body limit are not satisfied. As 
a result, the interaction of such small-size configurations is not tamed. Thus, in contrast to 
the Gribov approach, only a fraction of all configurations will interact according to the BBL 
approximation and therefore C < 1 in ([13]). 

Using ([!]), it is straightforward to estimate the kinematical boundaries where the unitarity 
limits may be reached. Indeed, the requirement that a e \ < atot/2 [0,0 (also assuming 
that (P is applicable for the range of x in which the taming is necessary) indicates that for 
some gluon-induced hard processes the unitarity limit should be well within the reach of an 



electron- nucleus collider at HERA/THERA (for a review see ||11|| ). 

One can also make an interesting prediction about nuclear shadowing. Since the in- 
dependence of the nuclear structure function of (|TT|) is significantly weaker than that for the 
proton structure function of fllPl) , nuclear shadowing is not saturated (as is often assumed 
in the limit of fixed Q 2 and x — > 0), and we find 



F 2 A (x,Q 2 ) _ R 



AF 2 N (x, Q 2 ) * AR 2 N 1 + Ac 2 N R\ In 2 l/x ' (15) 

However, at unrealistically small x, where the impact parameters become significantly larger 
than Ra, this effect will disappear. 
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In summary, since the contributions of small configurations (which have not reached 
the black body limit) remain significant in a wide range of x and Q 2 , studies of the total 
cross sections are a rather ineffective way to the search for the onset of the BBL regime. In 
particular, it may be rather difficult to distinguish the BBL from DGLAP approximation 
with different initial conditions in this way. 

In what follows we shall demonstrate that studies of DIS final states provide a number of 
clear signatures of the onset of the BBL regime, which will be qualitatively different from the 
leading twist regime. For simplicity, we will assume that the BBL is reached for a significant 
part of the cross section and, hence, restrict our discussion to DIS on a large nucleus so that 
edge effects (which are important in the case of scattering off a nucleon) can be neglected. 



B. Diffractive final states 

The use of Gribov's orthogonality argument (to neglect non-diagonal transitions in flT5|)) 
allows the integrals over the masses to be removed in the expressions for the structure 
functions and, since diffraction is 50% of the total cross section in the BBL, we immediately 
find for the spectrum of diffractive masses: 

dF° (3 \x,Q 2 ,M 2 ) ixR 2 A Q 2 M 2 p(M 2 ) 



dM 2 ~ 12tt 3 (M 2 + Q 2 ) 2 ' 

dF° {3 \x,Q 2 ,M 2 ) ixR\ Q 4 p(M 2 ) 



(16) 



dM 2 12tt 3 (M 2 + Q 2 ) 2 ' 

Moreover, the spectrum of hadrons in the centre of mass of the diffractively produced system 
should be the same as in e + e~ annihilation. Hence, the dominant diffractively-produced final 
state will have two jets with a distribution over the centre of mass emission angle proportional 
to 1 + cos 2 9 for the transverse case and sin 2 9 for the longitudinal case: 

dF^ 3 \x,Q 2 ,M 2 ) 3 , , p ^R 2 A Q 2 M 2 p{M 2 ) 

dM 2 dcos9 8 l J 12vr 3 (M 2 + Q 2 ) 2 ' 1 ' 



dF^ 3 \x,Q 2 ,M 2 ) _3 , rrR 2 A Q 4 p(M 2 ) 

dM 2 dcos9 4 12-7T 3 (M 2 + Q 2 ) 2 1 ' 
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The transverse momentum of the jet, pt, and the longitudinal fraction of photon energy, 
z, carried by the jet are related to the diffractive mass, M, and the angle, 9, as follows (we 
neglect here the quark masses as compared to Q, M): 

M ■ a 

p t = — - sin , 

z = (1 + cos#)/2 . (19) 
Hence, in the BBL, diffractive production of high p t jets is strongly enhanced: 

(p 2 t (jet)) T = 3M 2 /20, 

(p 2 tUet)) L = M 2 /5 . (20) 

This is to be compared to the leading twist approximation where it is oc InQ 2 . The relative 
rate and distribution of the jet variables for the three jet events (originating from qqg 
configurations) will be also the same as in e + e~ annihilation and hence is given by the 
standard expressions for the process e + e~ — > qqg (see e.g. 



An important advantage of the diffractive BBL signal is that these features of the diffrac- 
tive final state should hold for M 2 < Q\bl even for Q 2 > Q\bl because configurations with 
transverse momenta < Qbbl/^ still interact in the black regime (and correspond to trans- 
verse size fluctuations for which the interaction is already black). 

Another interesting feature of the BBL is the spectrum of the leading hadrons in the 
virtual photon fragmentation region produced in DIS. The spectrum is essentially given by 
the ^-dependence of (|T7| ) and ([18]). For fixed M 2 , the jet distribution over z for transversely 
polarised photons is simply 

^ocl + (2z-e (21) 
Similarly, for longitudinally polarised photons, 

^oc z(l-z). (22) 

dz 

If no special separation procedure is undertaken, at small x one actually measures <jl + e^T 
(e is the photon polarisation). In this case, combing ([[7]) and ([18]) we find: 



Exclusive vector meson production in the BBL is in a sense a resurrection of the original 
vector meson dominance model |^2| without off-diagonal transitions. The amplitude for 
the vector meson-nucleus interaction is proportion to 2tcR\ (since each configuration in the 
virtual photon interacts with the same BBL cross section). This is markedly different from 



the requirements [43 for matching generalised vector dominance model (see e.g., |Q) with 
QCD in the scaling limit, where the non-diagonal matrix elements are large and lead to 
strong cancellations. Hence, we can factorize out the cross section for the dipole interaction 
from the overlap integral between wavefunctions of virtual photon and vector meson to 
obtain for the dominant electroproduction of vector mesons 



dt Q 2 dt 

(2nR 2 A ) 2 3T v M* 4 \j x (y/=tR A ) 



2 



(24) 



16tt a(M 2 +Q 2 ) 2 -tR\ 

where T v is the electronic decay width V — > e + e~, a is the fine-structure constant, and 
J\{x) is the Bessel function. Thus the parameter-free prediction is that in the BBL at large 
Q 2 vector meson production cross sections have a 1/Q 2 behaviour, in stark contrast to an 
asymptotic behaviour of 1/Q 6 predicted in perturbative QCD |7|,p^,pJ[| . 

In order to observe the onset of the BBL regime for a nucleon target, one should consider 
scattering at small impact parameters. However, a direct comparison of the BBL prediction 
for small b with data is very difficult since, in the BBL, the amplitude oscillates as a function 
of t. Also, in the case of a proton target, a significant nucleon spin-flip amplitude may mask 
these oscillations. 



C. Inclusive spectra 



In the leading twist approximation, the QCD factorization theorem is valid and leads to 
universal spectra of leading particles (independent of the target) for the scattering off partons 
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of the same flavour. Fundamentally, this can be explained by the fact that, in the Breit frame, 
the fast parton which is hit by the photon carries practically all of the photon's light-cone 
momentum (z — > 1). Due to QCD evolution, this parton acquires virtuality ~ Q 2 , and a 
rather large transverse momentum, kt (which is still -C Q 2 )- So, in pQCD, quark and gluons 
emitted in the process of QCD evolution and in the fragmentation of highly virtual partons 
together still carry all the photon momentum. In contrast, in the BBL, configurations 
in which partons carry all of the photon momentum form only part of cross section of 
leading hadron production. Another part is from inelastic collisions of configurations where 
all partons carry appreciable momentum fractions and large relative transverse momenta 
(these configurations are rather similar to the case of diffractive scattering (see e.g. (|17D 
and flT8|))). Hence, in the BBL case, the spectrum of leading hadrons (in the direction of 
the virtual photon) is expected to be much depleted. 

The inclusive spectrum of leading hadrons can be estimated as due to the independent 
fragmentation of quark and antiquark of virtualities > Q 2 , with z and p± distributions 
given by ( |TTD and ( |T8D (cf. the case of diffractive production of jets discussed above). 
The independence of fragmentation is justified because large transverse momenta of quarks 
dominate in the photon wave function (cf. eqs. ( |TT| , [TS|, (TUJ) and because of the weakness of 
the final state interaction between q and q, since the a s is small and the rapidity interval is of 
the order of one. Obviously, this leads to a gross depletion of the leading hadron spectrum as 
compared to the leading twist approximation situation where leading hadrons are produced 
in the fragmentation region of the parton which carries essentially all momentum of the 
virtual photon f\ Taking the production of multi-jet states like qqg into account will further 
enhance this scaling violation. If we neglect gluon emissions in the photon wave function, 



Qualitatively, this pattern is similar to the one expected in the soft region since the spec- 
trum of hadrons produced in hadron-nucleus interactions is much softer than for hadron-nucleon 
interactions. 
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we find, for instance: 



dz 



J D^ h (z/y,Q 2 )±(l + (2y-l) 2 )dy, (25) 



for the total differential multiplicity of leading hadrons produced by transversely polarised 
virtual photons, in the BBL. D q ^ h (z/y,Q 2 ) is the fragmentation function of a quark with 
flavour q into the hadron h. Here we use that D u ^ h (z/y, Q 2 ) = D d ^ h (z/y, Q 2 ) and neglect a 
small difference in the fragmentation functions of light and heavy quarks. 

An illustration of the results of the calculation of dN^I 11 j dz is presented in Fig. |l|. We 
normalise the distribution to the leading twist case by using realistic up-quark fragmenta- 
tion functions at Q 2 = 2 GeV 2 |45| (the up and down quark distributions are similar and 
we neglect the small difference induced by including further flavours). One can see from 
the figure that a gross violation of the factorization theorem is expected in the BBL. The 
spectrum of leading hadrons is much softer at large z, with an excess multiplicity at z < 0.1. 
Note that the use of leading twist fragmentation functions in the above expression probably 
underestimates absorption. So the curve in Fig.|I] can be considered as a conservative lower 
limit for the amount of suppression. 

With an increase of Q 2 we expect a further softening related to a change in the partonic 
structure of the virtual photon wave function. Progressively more configurations contain 
extra hard gluons, which fragment independently in the BBL, further amplifying deviations 
from the standard leading twist predictions. 

Another important signature of the BBL is the change of pt distributions with decreasing 
x (at fixed Q 2 ). The spectrum of the leading hadrons should broaden due to increased p t 
of the fragmenting partons. Hence, the most efficient strategy would be to select leading 
jets in the current fragmentation region and examine the z and p t -dependence of such jets. 
Qualitatively, the effect of broadening of pt-distributions is similar to the increase of the p t 



distribution in the model ||46|| , although final states in DIS were not discussed in this model. 

An important advantage of inclusive scattering off a nucleus is the possibility to use 
a centrality trigger. For example, one could use the number of nucleons emitted in the 
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nucleus decay (soft nucleons in the nucleus rest frame) to select scattering at the central 
impact parameters. Such a selection allows the effective thickness of the nucleus to be 
increased, as compared to the inclusive situation, by a factor ~ 1.5 and, hence, allows for 
the BBL to be reached at significantly larger x. The signal for the BBL will be a change 
of the spectrum with centrality of the collisions, in contrast to the LT case where no such 
correlation is expected. Note that the lack of absorption of leading particles in DIS off nuclei 



at the fixed target energies and Q > 2 GeV is well established experimentally, see e.g. (47 



To summarise, predictions for a number of simple final state observables in the black 
body limit are distinctly different from those made in the leading twist approximation and, 
hence, will provide model-independent tests of the onset of the BBL. 



IV. UNITARITY CONSTRAINTS FOR ELECTRON-NUCLEUS DIS 

While predictions for inclusive DIS and DIS final state observables in the BBL are dis- 
tinct, one still would like to determine the kinematical region where the regime of the BBL 
sets in. One way to address this issue it to examine the unitarity of the S'-matrix for the 
interactions of purely perturbative QCD fluctuations. 

A. The interaction of small colour dipoles with hadrons 

The dipole picture of electron-target DIS is valid when the lifetime of the fluctuation 
is much longer than the interaction time with the target. Equally, this relation may be 
expressed in terms of the coherence length, l co h, of the fluctuation (relative to the target 
radius). This length is given by the average longitudinal distances (Ioffe distances) in the 
correlator of the electromagnetic currents which determine the structure function F 2 (x, Q 2 ). 
A simple analysis shows that l co h ~ l/(2m^x). At HERA l co h can reach values of 10 3 fm 
(for moderate Q 2 only). As Q 2 increases, the scaling violations lead to a gradual reduction 
of l co h at fixed x and ultimately to the dominance of longitudinal distances ~ (for a 
discussion see (ill). 
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In the target rest frame, when l co h is sufficiently large, the virtual photon can fluctuate 
into a variety of partonic configurations containing various numbers of partons and involv- 
ing different transverse sizes. The simplest is a qq dipole, which dominates at very small 
distances. Since under SU(3) C transformations the quark and the anti-quark transform in 
the fundamental representation (like 3 and 3) the name "colour triplet" is often used for the 
qq pair (although overall its colour is of course neutral). The gg-dipole of a small diameter 
interacts with a hadronic target with an interaction cross section given by ([!]). Another 
important photon fluctuation is the one consisting of a quark, anti-quark and gluon, with a 
relatively large transverse momentum between the quark and the anti-quark. Such configu- 
rations effectively transform in the adjoint representation (effectively 8, 8) and so are known 
as the qqg colour octet dipole. There are of course other qqg configurations, for example 
those in which the gluon and either the quark or anti-quark have a large relative transverse 
momentum (these merely correspond to the 0(a s ) corrections to the colour triplet dipole) 
and other more general configurations which do not correspond to dipoles at all. 

For the case of scattering of colour octet dipoles off a target, the corresponding cross 
section is enhanced a colour factor (given by the ratio of the Casimir operators of SU(3) C ), 
C F (8)/C F (3) = 9/4 



~ octet /j2 \_ 9 * triplet/ -a n 
°pQCm a _L' x ) — ^ °pQCD l a ±? x ) 

3tt 2 

= — dla s (Q 2 )xG T (x,Q 2 ). (26) 

Both (|) and © predict cross sections steeply rising with increasing energy (driven by the 
rise of the gluon density with decreasing x). Fitting the energy dependence in the form 7 
o" p qcd(s,Q 2 ) oc s n( - Q2 \ one finds 

n (Q 2 = 4 GeV 2 ) « 0.2 , 



7 Such fit is useful in practical applications. Perturbative QCD predicts the behaviour oc a + 
blnl/x + cm 2 l/x for the HERA energy range since radiation of < 1-2 hard gluons is possible in 
the multi Regge kinematics at HERA. 
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n (Q 2 = 40 GeV 2 ) » 0.4 . (27) 

One of the manifestations of the behaviour predicted by fl2"7|) is the Q 2 -behaviour of the 
measured exclusive vector meson production. The interaction cross sections of (|I|) and ( |26| ) 
may be thought of as a complimentary description of the physics described by the leading 
log QCD evolution equations at small x. However, an accurate determination of the relation 
between Q 2 and the transverse size d± in these equations requires a next-to-leading order 



QCD analysis in this framework which has not been done yet. Numerical studies ||i9| . |3T 
based on matching of the <ij_-space and Q-space expressions for ctl(x, Q 2 ), lead to A ~ 9 — 10 
for a sufficiently broad range of Q 2 and x. With this choice of A, a good description of 
the recent inclusive DIS electron-proton data was obtained (with a suitable extrapolation 
to large d± |3T|) without any further fitting. As already mentioned, A is a logarithmically- 
decreasing function of the dipole size d±. In particular, at large values of d±, dj_ > 0.3 
fm, where connection between d\ and Q 2 is rather sensitive to non-perturbative effects, one 
expects a decrease of A with increasing d\. On the other hand, it was found that variations 
in A do not significantly effect values of (XpQCD^l, x). Hence, our following estimates of 
the unitarity constraints, which are made using ([!]) and (p6|), are insensitive to a possible 
decrease of A at large d 2 ± . 

B. Unitarity constraints 

The rapid increase of the cross sections given in ([I]) and (|26|) with decreasing x cannot 
continue forever, otherwise the unitarity of the S'-matrix will be violated (see (^)). The 
unitarity boundary Im f(s, b, d 2 j_) = 1 can be expressed in terms of the pQCD dipole cross 
section as 

^qcd = ^ = cr tot /2 , (28) 

and is applicable to any hadronic target. For a nucleus, with the atomic number A, o"pQ°c D 
cannot exceed its geometric limit ttR 2 ^. Thus, generalising ([!]) and (26) for a nuclear target, 
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we obtain the following kinematical restrictions imposed by the unitarity of the S'-matrix 
for x <C l/iRAm^: 



7T~ 



af(d 2 ± ,x) = '-dl [xG A (x,Q 2 )\ a s (Q 2 ) <ttR 2 a , 



aT(dl,x) 



3vr 2 



xG A (x,Q 2 ) a s {Q 2 ) <nR 



A > 



(29) 



where xG (x, Q 2 ) is the nuclear gluon density The kinematical boundaries following from 



these equations are presented as curves in the x-Q plane in Figs. 14-17 of JIT 



The unitarity constraints are even more stringent for DIS on nuclei at central impact 
parameters b, b < R A . This is essentially equivalent to scattering off a cylinder of the length 
2R A oriented along the reaction axis. Obviously, in this case, edge effects are suppressed 
and one also gains an additional factor of ~ 1.5 on the left hand side of (^) due to the 
increased density of nucleons. The nuclear gluon distribution at a given impact parameter 
b is introduced as j|9 



xG A (x, Q 2 , b)=A xG N (x, Q 2 ) f A (x, Q 2 , b) T A {b) , 



(30) 



where the function f A (x,Q 2 ,b) describes the amount of nuclear shadowing, T A {b) = 
J^oo dzp A (b, z) and / d 2 bT A {b) = 1. Note that nuclear shadowing is larger at central im- 
pact parameters then in the situation when one averages over all impact parameters. Now, 
the unitarity constraints for DIS on nuclei at central impact parameters immediately follow 
from (|29|) for colour triplet 



7T - 



1.5 x — r" 
3 



xG A (x,Q 2 ,b = 0) a s {Q 2 ) <nR 2 A 



(31) 



and colour octet dipoles 



1.5 x 



3tt 2 



xG A (x,Q 2 ,b = 0) a s (Q 2 ) <tcR 2 a . 



(32) 



The kinematical regions prohibited by these unitarity constraints (defined by x < x\i m , Q < 
Q e s) lie to the left of the curves in Figs. ^ and [3] (on the boundary Q = Q c q, x = xji m ) . For 



the nucleon gluon density we used the CTEQ4L pi parameterization evaluated at the scale 
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Q 2 = 4 GeV 2 . For each nucleus, we present scenarios with the highest and lowest shadowing 
(see Sect. |V|). The curves with more shadowing lie below the ones with less shadowing, 
for all Q 2 . To illustrate the trends given by fl31|), the curves are extended to the region 
x > l/(4i?AWjv) where, strictly speaking, the unitarity constraints (of ([H]) and fl32|)) should 
not be directly applied. 

As one can see from Figs. ||] and || effects associated with the unitarity constraints are 
expected, in a wide range of x and Q 2 , to be covered by THERA. Regardless of the nature 
of such effects, strong modifications of the gluon field in heavy nuclei (as compared to the 
incoherent sum of the nucleon fields) appear to be unavoidable. 

In summary, in the search for the BBL we observe that employing nuclear targets allow 
us to gain substantially in the region where higher twist effects become important. However, 
the presence of the leading twist nuclear shadowing reduces the magnitude of this gain. It is 
very important that in a wide range Q 2 the unitarity limit is reached at relatively large x so 
that lnQ 2 /AQ CD is comparable to hix^/x (where Xq ~ 0.05 is the starting point for the gluon 
emission in the lnl/x evolution). Hence, the diffusion to the small transverse momenta is 
likely to be a correction. Therefore non-perturbative QCD, with a large coupling constant, 
is unlikely to be relevant for the taming of the structure functions. 



V. NUCLEAR SHADOWING AND DIFFRACTION 



A long time ago Gribov [51J established an unambiguous connection between the cross 
section of small-t diffraction of a hadron off a nucleon and the amount of shadowing in the 
interaction of the same hadron with a nucleus, for the limit in which only two nucleons of the 
nucleus are involved. Applying Gribov's formulae to describe photon-deuteron scattering, 
the effect of nuclear shadowing for the total cross section can be expressed as 



8 The contribution of the real part of the diffr active scattering amplitude was neglected in [51] 
since it was assumed that the total cross section is energy-independent. 
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<r to t(eD)-2<rtot(ep) _ (l-A 2 )^f^U 1 
shad a(ep) (1 + A 2 ) (T-tot(ep) 8*7$ ' 



where A is the ratio of real to imaginary parts of the amplitude for diffractive DIS, and 
Rd is the deuteron radius. For small x, A may be as large as 0.2 4- 0.3, which results in 
(1 — A 2 )/(l + A 2 ) ~ 0.8 4- 0.9. For simplicity, we have neglected the longitudinal nuclear 
form factor (which leads to a cutoff of the integral over the produced masses, cf. e.g., ([13])) 



It is worth emphasising that (|33| ) does not require the dominance of the leading twist in 
diffraction. The only assumption, which is due to a small binding energy of the deuteron 
and which is also known to work very well in calculations of hadron-nucleus total and elastic 
cross sections, is that the nucleus can be described as a multi-nucleon system 0. Under these 
natural assumptions, one is essentially not sensitive to any details of the nuclear structure, 
such as short-range nucleon correlations, etc. 

A. Inclusive diffraction at HERA and predictions for nuclear parton densities 

The Gribov approximation has been applied to the description of shadowing in nuclear 
DIS for a long time. The first calculations were performed in |5^ , |55[] before the advent of 



HERA and, hence, required modelling of diffraction in DIS. This modelling was based on the 
QCD extension of the Bjorken aligned jet model and produced a reasonable description 
of the NMC data ||56||. More recently, an explicit use of the HERA diffractive data allowed 



an essentially parameter-free description of these data []57j to be provided. 

Two important features of the HERA inclusive diffractive data [|5^,|)IJ are the approxi- 
mate Bjorken scaling for the diffractive parton densities and a weak dependence of the total 



9 The condition that the matrix element <^4|T[J^(y), Jj,(0)]|j4 > (which defines the nuclear struc- 
ture function and total cross section) involves only nucleonic initial and final states is not so 
obvious in the infinite momentum frame. However, it is implemented in most of the light-cone 
models [35,53j. 
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probability of diffraction -Pdiff 011 Q 2 '- -Pdiff ~ 10% at Q 2 > 4 GeV 2 . The first observation 
is in line with the Collins factorization theorem |60[ which states that in the Bjorken limit, 
the diffractive structure functions fP((3, Q 2 , Xp, t) satisfy the DGLAP evolution equations 
{(3 — x/xp « Q 2 /{Q 2 + M 2 ) and xp « (M 2 + Q 2 )/2q ■ (p — p') is the fraction of the proton's 
momentum carried by the diffractive exchange). 

A relatively small value of the probability of diffraction (as compared to the case of ttN 
scattering) indicates that the average strength of the interaction leading to diffraction is 
correspondingly small. To characterise this strength it is instructive to treat diffraction in 



the ^-channel picture using the eigenstates of the scattering matrix ||61|| . Such an approach 
is complementary to the picture of the factorization theorem. The use of the optical theorem 
leads to the following definition of the effective strength of the interaction (for diffraction in 
nucleon collisions and shadowing in nuclear collisions): 

^diffi 

aUx, Q 2 ) = 16tt . dt k=0 . (34) 

Here the superscript j indicates which hard parton is involved in the elementary hard process 
of the cross section atot(x,Q 2 ). This equation allows the average cross section for configu- 
rations which contribute to quark-induced and gluon-induced diffraction to be extracted. 

Recently a number of phenomenological analyses of HERA inclusive diffractive and 
diffractive jet production data were performed within the leading twist approximation for 
hard diffraction. We have compared several of these analyses [p8| , |62| , |63| and found that they 
lead to fairly similar values of a 3 ^, especially if a matching at large diffractive masses to the 
soft factorization was implemented (for details see ||64|). 

The results for a 3 eS for up quarks and gluons are presented in Fig. [|. We used results 



of two analyses of the HERA diffractive data (by the HI collaboration |58 ], labelled by 



"HI" in the figure, and by Alvero et al. |p2 |, labelled by "ACWT+"). The diffractive parton 
densities of |)2| were extended to the region of small f3 by requiring consistency with the soft 
factorization theorem. Three curves (solid, dashed and dotted) represent the Q 2 -evolution 
of a 3 e Q. While the Q 2 -evolution is not very significant for <jg ff of the up quark, it decreases 
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<7gfj rather rapidly for gluons. This rapid change is explained by large scaling violations for 

xG N (x,Q 2 ). 

Combining the Gribov theory with the Collins factorization theorem and comparing the 
QCD diagrams for hard diffraction and for nuclear shadowing arising from the scattering 
off two nucleons (see Fig. |5|), one can prove that, in the low thickness limit, the leading 
twist nuclear shadowing is unambiguously expressed through the diffractive parton densities 
f^(x/x]p, Q 2 , Xjp, t) of ep scattering: 

\ \ \2 p pOO pOO pXQ 

fj/A(x, Q 2 )/A = fj/ N (x, Q 2 ) - - 2 / d 2 b / dzi / dz 2 / dx F 

Z 1 T A J J—oo J z\ Jx 

xp A (b, z^paQ), z 2 ) cos(x 1P m N (z 1 - z 2 ))ff /N [p, Q 2 , x F , tj \ k 2 =0 , (35) 

where fj/A(%, Q 2 ) and fj/^(x, Q 2 ) are the inclusive parton densities, pa is the nucleon density 
in the nucleus. At moderately small values of x, one should also add a term related to 
the longitudinal distances comparable to the inter-nucleon distances in the nucleus. This 
additional term can be evaluated using information on the enhancement of the gluon and 
valence quark parton densities at i ~ 0.1 at the initial scale Qq. This would slightly 
diminish nuclear shadowing at higher Q 2 via the Q 2 evolution. To summarise the results of 
(P5j), the nuclear shadowing effect given by (1 — fj/A(x,Q 2 )/Afj/iy(x,Q 2 )), is proportional 
to a J eS (x,Q 2 ). 

For N > 3 nucleons, we need to establish which configurations in the photon wavefunction 
dominate in diffraction. There are two extreme alternatives, i.e. a dominance of either 
hadronic-size configurations (as indicated by the aligned jet model) or of small size (~ 
1/Q) fluctuations. Figure |] represents the corresponding diagrams (contributing to nuclear 
shadowing) . 

One can determine which of these extremes is closer to reality by comparing the effective 
cross section characterising diffraction ([34]) with the cross section, <T p Q CD ((i^, x), for the 
double interaction of a dipole with a small and fixed diameter d^ ^ 1/Q, within the eikonal 
approach, given by (|I|). Thus, for the x-range studied at HERA at Q 2 > 4 GeV 2 (x ~ 0.001) 
, large-size configurations dominate, while the small dipoles contribute little to the bulk of 
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the inclusive diffractive cross section. This can also be seen from a comparison of the amount 
of nuclear shadowing in calculated within the leading twist and eikonal approaches. This 
is consistent with the experimental observation that the t-slope of the inclusive diffraction 
[B ~ 7 GeV~ 2 ) is significantly larger than for the processes where small-size dipole dominate 
(B ~ 4.5 GeV -2 ). At lower x ~ 10 -4 , small dipoles may become much more important due 
to the increase of <7 p qcd due to the QCD evolution. 

It is well known that total, elastic and inelastic diffractive cross sections for the inter- 
actions of hadrons with nuclei are quantitatively well- described within the scattering eigen- 
state approximation [[n|], which is a generalization of the eikonal approximation. Since, at 
Q 2 = 4 GeV 2 , relatively soft interactions give the dominant contribution to diffraction, in- 
teractions with N > 3 nucleons can be considered using a generalised eikonal model with an 
effective cross section given by ([34]). We found that, due to a relatively small value of crj ff , 
the amount of nuclear shadowing for is relatively modest (as compared to the nuclear 
gluon distribution) and is practically independent of fluctuations of the effective cross section 
(when cr e ff is kept fixed) pEfl . The Q 2 -dependence of shadowing is taken into account via the 
DGLAP evolution equations. Recently, we performed an extensive comparison of different 
parametrisations of the diffractive structure functions p3 |. The results will update the anal- 



ysis |)5| to include a range of modern parameterizations of the quark and gluon diffractive 



parton densities. Using the effective cross sections presented in Fig. |J and generalising fl35|) 
to include the rescattering terms with N > 3 nucleons, we are able to produce predictions 
for up quark and gluon parton distributions for a number of nuclei. As an example Fig. [7] 
represents the ratio of nuclear to nucleon up quark and gluon parton distributions for nuclei 
of 12 C (carbon) and 206 Pb (lead). Note that since the models we used provide good fits to 
the HERA diffractive data for F®(x, Q 2 ,xjp), they effectively take into account higher twist 
terms, if they are present in the data. 

Hence, we conclude that combining the Gribov theory with the factorization theorem for 
inclusive diffraction and experimental information from HERA provides reliable predictions 
for nuclear parton distributions and, hence, for F£ in the HERA and THERA experimental 
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ranges. 



B. Gluon shadowing 

We explained above that since the interaction of the colour octet dipole is enhanced by a 
factor 9/4 relative to a colour triplet P2,25,4"8|, one expects, for processes sensitive to such 



configurations, an earlier onset of the regime where unitarity effects may become important. 

The leading twist mechanism of nuclear shadowing is connected to the amount of the 
leading twist diffraction in gluon-induced reactions. The larger strength of the perturba- 
tive interaction as well as stronger non-perturbative interactions up to a scale ~ 2 GeV 
in the gluon channel suggest that the gluon-induced diffraction should occur with a large 
probability. 

This is consistent with another interesting feature of the current HERA diffractive data 
(shared by the diffractive data from the proton colliders), i.e. a very important role of the 
gluons in hard diffraction. In the language of the diffractive community, the "perturbative 
Pomeron" is predominantly built of gluons. To quantify this statement it is convenient to 
define the probability of diffraction for a hard probe which couples to a parton j [ 65| : 



rj / ^ S 'dtdxp fP(x/x P ,Q 2 ,x P ,t) 

P > (X ' Q) - M ' (36) 

A large probability of diffraction in the gluon channel: 

P g (x < 10~ 3 ,Q 2 = 4 GeV 2 ) ~ 0.4, 
P g (x < 1(T 3 ,Q 2 = 10 GeV 2 ) ~ 0.2, 

(37) 

leads via fl36|) to a large (see Fig. (4|) and, hence, to the prediction of a large leading 
twist shadowing for G A (x < 10 _2 ,<5 2 ) (see Fig. 0). Interactions of small dipoles are more 
important in the gluon case than in the quark case. Nevertheless, if we use the eikonal 
approximation (in the same spirit as for qq dipoles) to estimate the relative importance of 
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the leading and higher twist effects, we find that up to fairly small x in a wide range of 
Q 2 , the leading twist terms dominate. Note that in this case one should use cr q qN from |JT 



rescaled by the factor 9/4 even for the dipole sizes where non-perturbative effects could be 
important. A comparison of Figs. |2] and ^ suggests that the higher twist effects appear to 
be more important in the gluon channel than in the quark channel. Also, the gluon induced 
interactions of a projectile with several nucleons are much more important and more sensitive 
to details of the interaction dynamics. However, with (j 9 ^ ~ 40 mb, a large leading twist 
shadowing in the gluon channel appears to be unavoidable. It seems to be large enough to 
reduce strongly the gluon densities. However, even this strong reduction can not prevent the 
violation of unitarity for the interaction of colour octet systems with heavy nuclei, as can 
be seen from Fig. |3|. For scattering at central impact parameters, this may already occur at 
Q ~ 2 GeV for the whole x range, x < 1/ {Aiti^Ra)- 



C. Shadowing in the interaction of small dipoles with nuclei 

For several small-x processes, we can probe the interactions of a small colour dipole with 
the nucleus directly. These include <tl(x, Q 2 ), (7(7 + A — > J/ip + A) and c*{j* L + A — > p + A). 

We argued above that the eikonal approximation for fixed small d± leads to much weaker 
absorption effects than the leading twist shadowing. The eikonal approach is also inconsis- 
tent with the QCD factorization theorem for the production of vector mesons |§ which leads 
at small x to R 



%hiA^VA)\^ 
(Y L N^VN)\ t=Q 



da 
dt 



Ffa, Q) 

F N< 



1 2 



xG A (x,Q) 
xG N {x,Q) 



(38) 



Another observable sensitive to the difference between the leading twist and eikonal 
approaches is F£. It can be expressed via nuclear parton distributions (we consider the 
leading order in a s expression): 

Fi(x,Q 2 ) = ^^J 1 ^(x/y) 2 (^e 2 (l-x/y)yG A (y,Q 2 ) + lF A (y,Q 2 )) . (39) 

^ x y i=i 

The amount of nuclear shadowing for o A can be expressed by the ratio 
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F A Si ^/y)Ha - x/y) J2i el yG A (y, Q 2 ) + 2F 2 A (y, Q 



AF l AJ^(x/ y y((l-x/y)Eiel yG N {y, Q 2 ) + 2F£ f (y, Q 2 )/3 



(40) 



Results using ( |4"0"D are presented in Fig. [S] for the same parameterizations as in the previous 
figures. 

The impact-parameter eikonal approximation for x <C l/(4mjviiU) (where finite coher- 
ence length effects, which cannot be unambiguously treated in the eikonal approximation, 
can be safely neglected), leads to substantially smaller shadowing for the interaction of small 
dipoles (see Fig. (15) of 1151). The cleanest way to observe this effect would be to study coher- 
ent J/ip production. THERA would be able to cover a broad range of x — (M 2 ,^ + Q 2 )/W 2 
starting from the region of colour transparency (amplitude proportional to nucleon number 
A) and extending to the colour opacity regime where there is a strong reduction of the 
amplitude (see Fig. |). 

D. Total cross section of coherent inclusive diffraction 

There is a deep connection between shadowing and the phenomenon of diffractive scatter- 
ing off nuclei. The simplest way to investigate this connection is to apply the AGK cutting 



rules [S7|. Several processes contribute to diffraction on nuclei: (i) coherent diffraction in 
which the nucleus remains intact, (ii) break-up of the nucleus (without producing of hadrons 
in the nucleus fragmentation region), (iii) rapidity gap events (with hadron production in 



the nucleus fragmentation region). In [52| we found that for x < 3 ■ 10 3 ,Q 2 > 4 GeV 2 , 



the fraction of the DIS events with rapidity gaps reaches about 30-40% for heavy nuclei, 
with a fraction of the events of type (iii) decreasing rapidly with A. Recently, together with 
M. Zhalov one of us (MS) investigated the dependence of the fraction of the events due to 
coherent diffraction and due to the break-up of the nucleus on the strength of the interaction, 
a eS' ^ e f° un d that this fraction increases with a 3 a rather slowly. Thus, it is not sensitive to 
fluctuations of <j c q. One can see from Fig. |l^ that one expects a significantly smaller fraction 
for quark-induced processes (~ 35%) of coherent events than for gluon-induced processes 
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(~ 45%). We also found that the ratio of diffraction with and without nuclear break-up is 
small (10-20%) in a wide range of nuclei, and weakly depends on cr e g . Hence, it would be 
pretty straightforward to extract coherent diffraction by simply using anti-coincidence with 
a forward neutron detector, especially in the case of heavy nuclei (see discussion in [^8] ) . 

It is worth emphasising that the proximity of cr e j/crt t to 1/2 does not necessarily imply 
the proximity of the BBL in the sense of hard physics (it may also occur in the soft aligned-jet 
type scenario). A key distinction in this regard is in the dominance of the dijet production, 



broadening of the p t distributions, etc, as discussed in Sect. pi . 

Other manifestations of leading twist shadowing include a strong A-dependence of the 
fluctuations of the central multiplicity distribution of the produced hadrons |52| , and strong 



modifications of the leading hadron spectrum |H| for rapidities where the diffractive con 



tribution is small. This is also in marked contrast to the black body limit, where major 
modifications occur already at the highest rapidities. 

VI. EXPERIMENTAL CONSIDERATIONS 

The expected nuclear effects in the small-x region are large enough so that measure- 
ments of the ratios of the corresponding quantities for electron-nucleus and electron-nucleon 
scattering with an accuracy of only a few percent will be sufficient. Hence, the luminosities 



required for the first run of measurements are rather modest. The estimates of |69| indicate 
that luminosities of about 1 pb _1 per nucleus would be sufficient for measurements of a 
number of important inclusive observables which include ratios of the structure functions, 
the A-dependence of global features of diffractive cross sections, ratios of the leading particle 
spectra and fluctuations of the multiplicity in the central rapidity range. Luminosities of the 
order of 10 pb _1 per nucleus will be necessary to measure gluon densities via 2 + 1 jet events 
directly, to observe the pattern of colour opacity in exclusive vector meson production and 
to establish differential diffractive parton densities of nuclei. 

Requirements for the detector are practically the same as for the study of small-x 
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electron-proton scattering. The only additional requirement is for forward neutron detectors 
with specifications similar to that for the ZEUS forward neutron calorimeter (FNC). These 
will be needed to trigger on central inelastic collisions by selecting events where a larger 
than average number of neutrons is produced. Note that for small x collisions at central 
impact parameters the virtual photon interacts inelastically with ~ A 1 ^ 3 nucleons. On the 
basis of an analysis |58f of E-665 data on the production of soft neutrons at x ~ 0.03 |70 



we 



estimate that the average number of soft neutrons produced in the central collisions will be 
~ 20 (as compared to ~ 4 for peripheral inelastic collisions). Overall, a strong correlation 
between the number produced nucleons and centrality of the event is expected at small x 
(as is also the case for hadron- nucleus scattering). As we explained above, such detectors 
would allow the breakdown of leading twist QCD to be studied in runs with a single heavy 
nucleus. 

A FNC- type neutron calorimeter would have nearly 100% acceptance for the neutrons 
from nuclear break-up. Since the average number of such neutrons exceeds one for A > 16, 
a FNC-type detector would allow an effective tagging of nuclear break-up. Hence, it would 
also simplify studies of diffractive processes with nuclei. The separation of the three classes 
of diffractive events (coherent scattering, nuclear break-up producing nuclear fragments 
(mostly neutrons and protons) in the forward direction, and diffractive dissociation with 
meson production in the nucleus fragmentation region from non-diffractive reactions) is very 
similar to the selection of diffractive events in the ep case (see the previous subsection). The 
contribution of diffractive dissociation to the overall diffractive cross section is significantly 
smaller than in the ep case. The experimental signatures for this class of reactions are similar 
to ep scattering since the energy flow in the forward direction is expected to have almost 
the same topology. Calorimetric coverage down to 9 < 1° in the forward region will allow 
for the detection of most dissociative reactions. Nuclear break-up is expected to constitute 
about 10% of coherent diffraction for a wide range of model parameters (see Fig. |I~0| ). In the 
case of coherent exclusive processes, such as coherent vector meson production, detection 
of the process will be further simplified since the average transverse momenta of coherently 
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produced vector mesons are p t < 2/Ra> which is much smaller than for exclusive processes 
with protons. 

VII. CONCLUSIONS 

The nuclear program at THERA has a very strong potential for the discovery of a 
number of interesting new high energy (small x) phenomena, including quark and gluon 
shadowing (for kinematics in which longitudinal distances are much larger than the heavy 
nuclear radius), hard coherent diffraction off nuclei and colour opacity in the production 
and interaction of small colour singlets. Crucially, for the first time, it will provide several 
effective tools for unambiguously probing the black body limit of photon-nuclear collisions, 
thereby investigating QCD in a new regime of strong interactions with a small coupling 
constant. 

We thank A. H. Mueller for a useful discussion and GIF, ARC, PPARC and DOE for 
support. LF and MS thank INT for hospitality during the time this work was completed. 



36 



FIGURES 
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FIG. 2. Unitarity boundaries for the interaction of the colour triplet dipole with nuclei at 
central impact parameters. Regions to the left of each curve are prohibited by the unitarity bound 
of (31). Two sets of curves are given for each nucleus correspond to two different models of leading 
twist shadowing (as discussed in subsection |V|) . 
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FIG. 3. Unitarity boundaries for the interaction of colour octet dipole with various nuclei at 
central impact parameters. Regions to the left of each curve are prohibited by the unitarity bound 
of (|32|). Two sets of curves are given for each nucleus correspond to two different models of leading 
twist shadowing (as discussed in subsection |V|) . 
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FIG. 4. The effective cross section for nuclear shadowing (or nucleon diffraction) of (34), er^, 
for a parton of type j = u,g, as a function of Bjorken x at several values of Q 2 . The solid curve 
is for Q = 2 GeV, the dashed curve is for Q = 5 GeV, the dotted one is for Q = 10 GeV. Two 
representative [63] sets of diffractive parton distributions have been used in the curves labelled by 
"ACWT+" || and "HI" @ (see text). 



40 




FIG. 5. Diagrams demonstrating the relationship between hard electron-proton diffraction, 
involving a parton of type j, and nuclear shadowing in inclusive DIS. Here the Pomeron symbol, 
IP, merely represents a generic label for vacuum exchange. 
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FIG. 6. Typical diagrams for the leading twist (left) and eikonal (right) models for nuclear 
shadowing involving a small dipole and two nucleons. The curly lines represent gluons. In the 
diagram on the right the blob and its attached legs represent the gluon distribution of the nucleon. 
The Pomeron symbol, IP, in the diagram on the left merely represents a generic label for vacuum 
exchange. 
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FIG. 7. The ratio of nuclear to nucleoli up-quark and gluon parton distributions as a function 



of x, scaled by nucleon number, A. Two representative [33 sets of diffractive parton densities, 
"ACWT+" and "HI" |f)8| , are used to calculate nuclear shadowing (see text) for each nucleon 
and taken together they give an indication of the spread of theoretical predictions. The solid, 
dashed, and dotted curves are for Q = 2,5, 10 GeV, respectively. 
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FIG. 8. The ratio of total cross sections for longitudinally polarised photons scattering on 

nuclear and nucleon targets, af; /(Acrf^), as function of x for Carbon and Lead. Two representative 
H| sets of diffractive parton densities are used, "ACWT+" Q and "HI" |58fl , to calculate nuclear 
shadowing. The solid, dashed and dotted curves are for Q = 2,5, 10 GeV, respectively. 
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FIG. 9. The colour opacity effect for the ratio, R, of the coherent production of J ftp and T 
from Carbon and Lead, normalized to the value of this ratio at x = 0.02, calculated in the leading 



twist models of gluon shadowing [65] both with (dashed curves) and without (solid curves) taking 
the fluctuations of the cross section into account. 
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FIG. 10. Fractions of the total cross section, given by inelastic scattering (in), coherent (el) 
and incoherent (q — el) diffractive scattering cross sections as a function of the average interaction 
strength for quarks and gluons (o^g of (|34|), cf. Fig. ^|) for scattering off Lead. 
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